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Abstract 

The hypothesis about the relation between the observed alignment of spots in the x-ray 
film in cosmic ray emulsion experiments and the features of events in which jets prevail 
at super high energies is tested. Due to strong dynamical correlation between jet axis 
directions and that between momenta of jet particles (almost collinearity) , the evaluated 
degree of alignment is considerably larger than that for randomly selected chaoticly located 
spots in the x-ray film. It appears comparable with experimental data provided that 
the height of primary interaction, the collision energy and the total energy of selected 
clusters meet certain conditions. The Monte Carlo generator PYTHIA, which basically 
well describes jet events in hadron-hadron interactions, was used for the analysis. 
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1 Introduction 



The intriguing phenomenon, the strong collinearity of cores in emulsion experiments [1], closely 
related to coplanar scattering of secondary particles in the interaction, has been observed long 
time ago. So far there is no simple satisfactory explanation of these cosmic ray observations in 
spite of numerous attempts to find it (see, for instance, [2] and references therein). Among them, 
the jet-like mechanism [3] looks very attractive and gives the natural explanation of alignment 
of three spots along the straight line which results from momentum conservation in a simple 
parton picture of scattering. Besides, the strong momentum correlation of particles inside a jet 
and correlation between jet axes due to singularity of QCD matrix elements allow us to suggest 
the high degree of alignment for more than 3 spots. This has been already demonstrated for 
the four cores in [3] but using a simplified picture of hadronization. 

With increasing energy of colliding hadrons (nuclei) hard and semi-hard jets begin to play 
an important role due to growth of their production cross sections. Thus the jet activity is likely 
to be a feature of all events above certain threshold collision energy. One of the manifestations 
of this activity and the strict momentum ordering inside a hard enough jet can be the observed 
strong collinearity of spots in emulsion experiments. The main purpose of the present paper 
is just to trace this relation in detail. In Sect. 2 we formulate the problem on the whole. 
Section 3 describes the results of numerical simulation made under conditions close to emulsion 
experiments in the framework of PYTHIA [4] , and some discussion. A summary can be found 
in Sect. 4. 

2 Problem under consideration 

In the Pamir experiment the observed events (7-hadron families with the alignment) are pro- 
duced, mostly, by a proton with energy > 10 4 TeV interacting at a height of several hundred 
metres to several kilometres in the atmosphere above the chamber [2, 5]. The collision products 
are observed within a radial distance up to several centimetres in the emulsion where the spot 
separation is of the order of 1 mm. One can estimate the typical transverse momentum in the 
events under consideration using the ratio (see also (4)): 

pxh = rE, (1) 

where E is the energy deposition in the spot, r is its spacing in the x-ray film, h is the height 
of interaction. So px is of the order of 10 GeV for r = 15 mm, h = 1 km, E = 700 TeV. The 
particles with such transverse momenta can be typically initiated by a jet with pj? ~ 50 GeV 
or larger, since the most probable value of a fraction of jet energy carried by leading particles 
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is ~ 0.2 [4]. Such energetic jets are already enough collimated: their effective angular cone 
size 9 e g < f 5° due to the strict ordering of transverse and longitudinal particle momenta in the 
leading logarithm approximation of perturbative QCD [6]. QCD "teaches" us also that this 
effective size decreases with the growth of the "jet hardness" (transverse momentum) as 

6 cS ~ (InfpfVAqcD))- 1 , (2) 

where Aqcd is the dimensional QCD parameter. 

The main conjecture is that the particle distribution from the hard enough jets in the x-ray 
film plane can lead to the alignment of the emulsion spots due to the strong collimation of such 
particles and dynamical correlation between jet axis directions. Let us consider the kinematics 
in detail. For our analysis it is convenient to parametrize 4-momentum of each produced 
particle i under consideration with its transverse momentum pxi (relative to the collision axis 
z), azimuthal angle fa and rapidity rji in the center-of-mass system: 

[\JpTi + mf cosh 7ft, p Ti cosfa, p Ti smfa, ^p 2 Ti + mf sinh^]. (3) 

In this case the transformation from the center-of-mass system to the laboratory one reduces 
to a simple rapidity shift: Q = Vo + Vu where Q, r/ are the rapidities of particle i and the 
center-of-mass system correspondingly in the laboratory reference frame. If we neglect the 
further interactions of particles propagating through the atmosphere (this gives the maximum 
estimation of the alignment effect), then their position in the transverse (xy)-plane is easily 
calculated 

Ti = — h = . h , (4) 

v *i yPri + m i sinh(?7o + ?7j) 

where v ri and v Z i are the radial and longitudinal components of particle velocity respectively 

(Ei = \jp\i + mf cosh(r]o + Vi) 1S ^ ne particle energy in the laboratory frame). 

Since the size of the observation region is of the order of several centimetres, these radial 

distances must obey the following restriction: 

r min < r ii (5) 

Ti < r max . (6) 

We set r min = r rcs ~ 1 mm, r max ~ 15 mm. The restriction (5) simply means that spots 
are not mixed with the central one formed by the particles which fly close to the collision axis 
(mainly from the fragmentation region of an incident proton). The separation of spots in the 
x-ray film gives another restriction on the distance between particles 



It must be larger than 1 mm: 

dij > T r es j (8) 

in the opposite case the particles must be combined in one particle-cluster until there remain 
only particles and/or particle-clusters with the mutual distances larger than r rcs , each such 
particle-cluster being considered as a single particle with coordinates defined in the same way 
as center-of-mass coordinates of two bodies: 

r ij = (T i E i + T j E j )/(E i + E j ). (9) 

Then we select 2, 7 clusters/particles which are most energetic and obey the restrictions (5, 
6, 8) and calculate the alignment Ajv c using the conventional definition [2, 5]: 

£fey fc cos(2<k jfc ) 
XN < ~ N C (N C -1)(N C -2Y (10) 

and taking into account the central cluster, i.e. N c — 1 = 2, 7. Here 4>ijk is the angle between 
two vectors (r k — rj) and (r k — r : ) (for the central spot r = 0). This parameter characterizes the 
location of N c points just along the straight line and varies from — 1/(JV C — 1) to I. For instance, 
in the case of the symmetrical and close to most probable random configuration of three points 
in a plane (the equilateral triangle) A 3 = —0.5. The ultimate case of perfect alignment is 
\ Nc = 1 when all points lie exactly along the straight line, while for an isotropic distribution 
Ajy c < 0. The alignment degree P^ c is defined as a fraction of events with \ Nc > 0.8 [2] with 
the number of cores not less than N c . 

3 Numerical results and discussion 

If the hypothesis about the relation of alignment to the prevailing jet character of events at super 
high energies is valid, then this must manifest itself first of all in nucleon-nucleon collisions. 
Therefore, to be specific we consider a collision of two protons and fix a primary energy in 
the laboratory system E\ ah ~ 9.8 x 10 4 TeV, that is equivalent to y/s ~ 14 TeV — just the 
energy attainable at LHC (the rapidity shift being r/ ~ 9.55 after the transformation from 
the center-of-mass system to the laboratory one). To simulate a collision of two protons with 
such energies we use the Monte Carlo generator PYTHIA [4], which basically well describes jet 
events in hadron-hadron interactions and is tuned using the available experimental accelerator 
data. 

The results of numerical simulation which follows the consideration in the previous Section 
are presented in Fig. 1 (solid curve) with the parameters r m i n = r rcs = 1 mm, r max = 15 mm, 
h = 1000 m, which are close to the conditions of emulsion experiments, with the additional 
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restriction on the energy threshold of particle registration in the emulsion: > E thr = 4 
TeV 1 . The estimated alignment degree P Nc for N c cores is considerably larger than that for 
randomly selected chaoticly located spots in the x-ray film, but is still too small (by a factor 
of 3 — 4) to describe the experimental data [7] even taking into account their large errors. This 
can mean that the jet activity is not sufficient at such energies or the jet mechanism can not, 
in principle, give the large experimentally observable alignment. 

In order to try to answer this question let us consider the influence of the applied restrictions 
(5, 6) (the laboratory acceptance criterion) on the spectrum of particles selected to calculate 
the alignment. For particles with high enough transverse momenta pt% relative to their masses 
mi these conditions (5, 6) reduce, mainly, to the restriction on the available particle rapidities 
in the center-of-mass system: 

r m in < Ti =>- T]i < 7] max = ln(r /r min ) ~ 4.95, (11) 

Ti < r max r)i > i] min = ln(r /r max ) ~ 2.25, (12) 
since in this case ~ r /e Vi for i] + rji > 1, where 

r = 2h/e r >°. (13) 

Due to the kinematical restriction [6], 

e^V? S V~s, (14) 

a production of harder jets with larger rapidities becomes possible with the growth of y/s. The 
rapidity region (11, 12) just corresponds to the transition from soft to hard QCD physics, where 
the jet activity could manifest itself. 

Here one should note that ultrarelativistic particles (pxi rn^) are detected in the x-ray film 
from the restricted rapidity region (11, 12) which excludes such configurations as back-to-back 
hard jets with rapidities close to zero in the center-of-mass system. But just such configurations 
with scattering of hard partons at angles close to 90° in the considered hadronic center-of-mass 
system (which in this case practically coincides with the partonic center-of-mass system) can be 
expected to be responsible for the alignment phenomenon. The point is that leading particles 
from both these hard jets have quantitatively comparable energies in the laboratory frame 
together with the "strong memory" of scattering plane. Meanwhile leading particles from any 
other back-to-back hard jets with the relatively large modulo rapidities | ± ^ 6 c * m . s .| in the 
partonic center-of-mass system have essentially different laboratory energies due to Lorentz 
: Our calculations are practically insensitive to this threshold in the wide interval of its varying. 
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boost. And the energy distinction is mainly determined by the value of exp(2| ± f]^l ms \). In 
the latter case particles from a forward hard jet produce as a rule the most energetic clusters 
(apart from the central one) in the laboratory frame, even if the particles from a backward (in 
the partonic center-of-mass system) jet hit the detection region. However such most energetic 
clusters from one jet are less correlated with the primary scattering plane and therefore will 
not be much aligned as clusters from both hard jets. This argumentation is confirmed by our 
simulation. 

In this connection it is necessary to comment on the work [3] in which the jet hypothesis has 
been suggested for the first time for the explanation of alignment phenomenon. There the high 
degree of alignment has been demonstrated for the four cores only, using a simplified picture of 
fragmentation process. In fact, an axis distribution has been calculated in the partonic center- 
of-mass system in the first order of perturbative QCD theory at the partonic level, considering 
three partons in the final state only. Then the Lorentz transformation has been done in order to 
find their directions and localizations with respect to the central spot in the laboratory frame, 
considering each parton-jet as one long-living system (the fragmentation time is of the order 
of flight time) with some effective mass and aggregate group velocity. The velocity has been 
fixed (f3* ~ 0.7 in [3]) so as to be able to include the events giving the high degree of alignment 
and corresponding to the two final-state parton-jets in the backward hemisphere and one in the 
forward hemisphere in the partonic center-of-mass system. In our variables this means that a 
mass factor p?/ \Jpt + M 2 must be very small for a such massive system (M ^> p T ) so that 
it hits the detection region even with negative rapidities (i.e. corresponds to the backward 
hemisphere in the partonic center-of-mass system), if one uses the same boost parameters as 
for the hadronic center-of-mass system without taking into account the possible additional 
boost due to the distinction between these partonic and hadronic frames. Note that for real 
particles, e.g. 7r-mesons which mainly contribute to the multiplicity, this mass factor becomes 
significant for very small transverse momenta, px-w m n = 0.14 GeV, only. However, as 
our investigation shows, falling of appropriately correlated particles into the observation region 
is still not sufficient to obtain the high degree of alignment because of the energy selection 
procedure, if the total number of particles is large and they generate many distinctly separated 
spots. 

For completeness one should also mention that high transverse momentum jet production 
has a connection to the double-core configuration of cosmic-ray events as it has been pointed 
out in [8]. Under certain conditions a hard forward (in the partonic center-of-mass system) jet 
together with a central bunch gives two relatively far separated clusters with large energies. 
The detailed studies of double-core (or binocular) phenomena with estimations of event rates 
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and average lateral spread of the 7-family using a PQCD based Monte Carlo can be found, for 
instance, in [9, 10]. 

Ultrarelativistic particles from the central rapidity region in the hadronic center-of-mass 
system (as possible sources of appropriately correlated spots) can hit the observation region 
owing to the decrease of r only, i.e. the decrease of the height h of primary interaction or the 
increase of the rapidity r] of the center-of-mass system due to the growth of energy y/s, as it 
follows from (13). The energy growth seems preferable, if we intend to be closer to emulsion 
experiments and increase the jet activity. However this demands the extrapolation of PYTHIA 
parameters and their special tuning to the experimentally untested energy domain. Updating 
can be done appropriately after the LHC operation starts. Moreover at present this generator 
already uses the extrapolation of experimentally tested cross sections and structure functions to 
the LHC energy region y/s ~ 14 TeV in order to estimate the effects expected at such energies. 

For illustration we utilize the first "less dangerous" alternative — decrease the interaction 
height by a factor of 20 rather than increase the energy y/s by the same factor of 20 at the 
initial height so that particles from both hard jets (with back-to-back structure), hitting the 
registration region, come from some rapidity range near rji ~ including adjoint positive and 
negative values. In this case the alignment degree becomes strongly dependent on the minimum 
transverse momentum of hard process, p^ ard , which is a parameter of PYTHIA. At the height 
h = 1 km such dependence was not visible, although we might catch some marginal tendency 
of the alignment degree to grow with the increase of j9^ ard at that height. However without 
the restriction on p^ rd from below (minimum bias) the result coincides practically with one 
obtained earlier (solid curve in Fig. 1) that shows some general characteristics of jet structure 
of events. If > 3 TeV, particles from these hard jets together with particles flying close 
to z-axis (within the transverse radius < 1 mm) result in the alignment degree (dashed curve) 
comparable with the experimentally observed one [7]. 

Thus the jet-like mechanism can, in principle, attempt to explain the results of emulsion 
experiments. For such an explanation it is necessary (but not sufficient) that particles from 
both hard jets (with rapidities near rji ~ in the center-of-mass system) hit the observation 
region. This is possible at the relatively small height h = 50 m and y/s ~ 14 TeV; or at 
the height h = 1000 m, but the considerably higher energy y/s ~ 14 x 20 = 280 TeV; or at 
some reasonable and acceptable intermediate combination of h, y/s and r max which meets the 
following condition: 

r = 2h/e r >° = 2hm p /y/s~ < kr max , (15) 

where m p is the proton mass. A;~l/2<lis needed in order to have particles with rji < that 
hit the detection region (see (12)). We verified the decisive significance of condition (15) to 
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allow the observation of large degree of alignment and its dependence on the process hardness 
for the smaller energy y/s ~ 1.4 TeV (where the prediction of PYTHIA is quite adequate) 
and the height h = 5 m (in accordance with (15)) thereby confirming this peculiar kinematic 
"scaling" . 

At p^ ard = 3 TeV jets carry away about half of the energy of colliding protons in the center- 
of-mass system due to the relationship in a parton picture £ ~ 2pj^/ y/s, where £ is a fraction 
of proton energy carried by each interacting parton (quark or gluon). The striking feature of 
such configurations in the rr-ray film is approximate equality of energy deposition in the central 
and the rest most energetic clusters, that can be one of the physical guideline to select the 
events with very hard jets not only at the generator level (simulation). If we simply apply 
the additional threshold on the minimum energy of detected clusters needed in the alignment 
analysis, then we still obtain neither the desirable selection of jet hardness nor the increase of 
the alignment degree. The small variation of resolution parameter r res does not provide the 
desirable effect also. However, introduction of another threshold on the total energy of all 
(N c — 1) selected clusters E% T ~ -Eiab/2 (without taking into account the energy deposition in 
the central cluster around r = 0), 

N c -1 

£ E i > E * r > ( 16 ) 
i=i 

allows us to select the events with hard jets only in a "natural" physical way and to reduce the 
hypothesis to the really active mechanism. Figure 2 shows that the alignment degree increases 
with the growth of (the restriction on p^ ard is absent at all!), and it becomes large enough 
(dashed curve) and comparable with the experimentally observed one [7] above the threshold 
Ef r ~ 0.1.Ei a b ^ 10 PeV. Though one should note that our estimations give still too steep 
dependence on N c as one can see in Figs, lb, 2b from comparison of slopes of straight lines 
with the experimental behaviour. 

To give the reader a feeling for the various measures of alignment we present in Figs. 3 
and 4 the spatial distributions of most energetic clusters in the (rcy)-plane for a few generated 
events along with the corresponding values of Aat c . Some spots are hardly visible because of 
their small sizes which are proportional to the cluster energies (especially in the case A4 > 0.8) 
or because they are outside a square 10 x 10 mmxmm (but inside a circle r = 15 mm) as 
it sometimes happens in the case As > 0.8. Besides for A4 > 0.8 we can distinctly see three 
relatively large spots resulted from two hard jets and a central bunch. 

Here one should note that there was slightly other criterion for the selection of families for 
the analysis in the works of Pamir Collaboration: the families with the total energy of 7-quanta 
larger than a certain threshold and at least one hadron present were selected and analyzed. The 
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alignment becomes apparent considerably at J2^-y > 0-5 PeV (the families being produced, 
mostly, by a proton with energy > 10 PeV). Since the adequate comparison of our estimations 
with experimental data is impossible without a full simulation of particle propagation through 
the atmosphere, taking into account the energy distribution of primary cosmic particles, etc., 
then in order to demonstrate the possibility of appearance of high alignment degree due to 
the jet mechanism we restrict ourself to the simpler (but as concerns physics essentially close 
to experiment) criterion of selection over the total energy of all particles. These particles are 
mostly 7r-mesons, the neutrals among them being the main source of the detected 7-quanta. It 
is natural that the threshold on the total energy of all particles must be larger than the similar 
threshold on the total energy of 7-quanta at the same collision energy. For comparison we 
estimate also the alignment degree selecting only the most energetic 7-quanta with their total 
energy larger than certain threshold E^ ir (Fig. 5): 

JVc-l 

E E h > K r - ( 17 ) 
1=1 

The result is close to that obtained previously with the threshold imposed on the total energy 
of all particles. 

Besides for jet events 

Pv ' = const (18) 



with a high accuracy (see Figs, lb, 2b, 5b, which present the dependence of alignment degree 
on the number of considered cores at the different values of hardness parameter (lb) and 
threshold total energy (2b, 5b) in the logarithmic scale). This constant depends on Pr xA , E% r ', 
E^ r , decreasing with their growth, and could in principle be determined by the kernels of the 
Gribov — Lipatov — Altarelli — Parisi — Dokshitzer equations [6, 11, 12, 13, 14] which describe 
the process of radiation of quarks and gluons in the initial and final states. And, in fact, this 
process is implemented in the PYTHIA generator together with the subsequent hadronization 
of quarks and gluons. 

If nevertheless particles from the central rapidity region rji ~ and the jet-like mechanism 
are insufficient to describe the observed alignment and there is another mechanism of its ap- 
pearance at the energy y/s ~ 14 TeV and the height h ~ 1000 m (mostly used in emulsion 
experiment estimations), then in any case some sort of alignment should arise at LHC too in 
the rapidity region (11, 12). This region must be investigated more carefully on the purpose 
to study the azimuthal anisotropy of energy flux in accordance with the procedure applied in 
the emulsion and other experiments, i.e. one should analyze the energy deposition in the cells 
of rj x 0-space in the rapidity interval (11, 12) (the equivalent threshold minimum particle 



8 



energy being Ef^ s = E thr / cosh 770 ~ 2E thr /e V0 ~ 0.6 GeV in the center-of-mass system). 
Note that the absolute rapidity interval can be shifted: it is necessary only that the difference 
(^max — ^min) is equal to ~ 2.7 in accordance with the variation of radial distance by a factor of 
15 (Vmax/rmin = 15) due to the relationship r\ ~ r /e Vi (independently of r ). In other words, 
since we use particle momenta in the center-of-mass system, then future data should be treated 
in accordance with the algorithm described earlier in Sects. 2, 3 introducing the corresponding 
laboratory observables. 

4 Conclusions 

Our analysis shows that for pp-collision at a fixed height of primary interaction above the 
energy y^i, when the condition (15) is fulfilled — that is ultrarelativistic particles from the 
rapidity interval near rji ~ in the center-of-mass system fall into the observation region inside 
the radius r max in the laboratory frame due to the large Lorentz factor — the alignment of 
spots arises (this, in principle, explains the existence of the experimental energy threshold of 
this effect) and the alignment degree becomes strongly dependent on the process hardness. If 
the process hardness is close to maximum for the given energy y/s, the estimated degree of 
alignment is already comparable with the experimentally observed one. Introducing another 
additional threshold (the scale of which is determined by the energy of an incident proton) on 
the total energy of all (N c — 1) selected most energetic clusters (without taking into account 
the energy deposition in the central cluster) allows us to select the events with high hardness 
in a "natural" physical way and thereby support the jet-like hypothesis, which later on may be 
accepted (or refuted) in further investigations of, for instance, the energy cluster distribution 
and their particle composition with regard for interactions in the atmosphere, etc. 

Meanwhile we suggest the more careful investigation of the rapidity region (11, 12) at 
LHC in order to reveal the new still unknown mechanisms of alignment if they exist. For 
this purpose one should perform the analysis of energy deposition in calorimeters of CMS and 
ATLAS experiments in accordance with the procedure described in Sects. 2, 3 (i.e. calculating 
the appropriate observables in the laboratory frame). Such investigation can clarify the origin 
of the alignment, test the alternative hypotheses and give the new restrictions on the values of 
height and energy. 
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Figure 1: The alignment degree Pn c as a function of cluster number N c at h = 50 m and -^/s = 14 
TeV in linear (a) and logarithmic (b) scales. The solid curve is the result (coincident with one at 
h = 1000 m) without restriction on the minimum value of process hardness pj? rd , the dotted curve - 
at pj? rd = 300 GeV, the dashed curve — at pj? rd = 3 TeV. Points (o) with errors are experimental 
data from [71. 
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Figure 2: The alignment degree Pn c as a function of cluster number N c at h = 50 m and -^/s = 14 
TeV in linear (a) and logarithmic (b) scales. The solid curve is the result (coincident with one at 
h = 1000 m) without restriction on the total cluster energy E^ 1 , the dotted curve — at E^ r = 2 
PeV, the dashed curve — at = 10 PeV. Points (o) with errors are experimental data from [7]. 
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Figure 3: Samples of core distributions for simulated events with E^ r = 10 PeV and A4 > 0.8. The 
size of spots is proportional to their energy (except for the central spot which is not to scale). 



13 





Figure 4: Samples of core distributions for simulated events with E^ r = 10 PeV and \% > 0.8. The 
size of spots is proportional to their energy (except for the central spot which is not to scale). 
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Figure 5: The alignment degree Pjy c as a function of cluster number N c at h = 50 m and ^/s = 14 
TeV in linear (a) and logarithmic (b) scales. The solid curve is the result (coincident with one at 
h = 1000 m) without restriction on the total energy of 7-quanta £'* hr , the dotted curve — at E^ r = 1 
PeV, the dashed curve — at £'* hr = 5 PeV. Points (o) with errors are experimental data from [7]. 
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